Coherent molecular bound states of bosons and fermions near a Feshbach resonance 
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We analyze molecular bound states of atomic quantum gases near a Feshbach resonance. A simple, 
renormalizable field theoretic model is shown to have exact solutions in the two-body sector, whose 
binding energy agrees well with observed experimental results in both Bosonic and Fermionic cases. 
These solutions, which interpolate between BEC and BCS theories, also provide a more general 
variational ansatz for resonant superfluidity and related problems. 
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The coherent transformation of a cold atomic gas to 
molecules in the vicinity of a photoassociation |]| or Fes- 
hbach resonance has enabled a fascinating probe of 
quantum dynamical behavior in coupled atom-molecular 
systems, together with remarkably precise measurements 
of quantum binding energies. Recent Bosonic experi- 
ments have extended the available species to 133 Cs, 87 Rb, 
and 23 Na [2j- Experiments on ultracold degenerate Fermi 
gases of 40 K and 6 Li atoms have resulted in spectacu- 
lar demonstrations of molecular Bose-Einstein conden- 
sate (BEC) formation 0, [j| and of possible fermion su- 
perfiuid behavior in the BEC-BCS crossover region 

Since these are many-body systems, it is useful to try 
to develop the simplest possible field-theoretic model that 
can explain their behavior. An essential feature of any 
correct many-body treatment is that the basic theory 
must be able to reproduce the physics of the two-body 
interactions. In this paper, we combine previous analytic 
solutions of a coherently coupled field theory 0, 
with an exact renormalization of the coupling constants 
[Io|. in order to obtain analytic predictions for the two- 
body bound states. This gives a unified picture of an 
Feshbach resonance experiment and related studies 0, 

H, E3, E3, E3, EE E3, EE EE E3, EE provided 

a small number of observable parameters are available. 
The predictions will be compared with experimental data 
and with coupled-channel calculations. 

To quantitatively model these experiments, consider 
an effective Hamiltonian for the molecular field (^o) in 
the closed channel and the atomic fields (^1(2)) in the 
free-atom dissociation limit of the open channel: 



Hi — Hq 



d 3 x 



(1) 



with the commutation (+) or anti-commutation (-) re- 
lation [*i(x,t),^t(x',t)]± = SijS(x - x') for Bosonic 

or Fermionic field operators respectively. The free 
Hamiltonian Hq includes the usual kinetic energy terms 
and the potential energies (including internal energies) 
due to the trap potential hVi(x), while Uij is the atom- 
atom, atom-molecule, and molecule-molecule coupling 
due to s-wave scattering. The atomic and molecu- 
lar masses are mi 2 and mo = mi + m.2, and E m 



11 [Vq(0) — Vi(0) — Va(0)] gives the bare energy detuning 
of the molecular state with respect to free atoms. 

Next, we consider a coherent process of Raman pho- 
toassociation or a magnetic Feshbach resonance coupling, 
giving rise to an overall effective Hamiltonian term in the 
homonuclear case (only with bosons) 0, IE 



H = Hi 



(2) 



or, for the case of heteronuclear dimer formation involv- 
ing either fermions or bosons 0]: 



H = Hi+H X 



dh 



(3) 



Here, x 1S the bare atom-molecule coupling responsible 
for the conversion of free atom pairs into molecules and 
vice versa. The heteronuclear case can be applied to 
Fermionic atom pairs in different spin states ($1, ^2) 
combining into a Bosonic molecule (^o), or pahs of 
Bosonic and Fermionic atoms combining into a Fermionic 
molecule, or else to a fully Bosonic case where the atom 
pairs are not identical. 

Bosonic homonuclear case. First we consider the fully 
Bosonic uniform case of Eq. J2J, i.e., a single-species 
atomic BEC (with mi = m) coupled to a molecular BEC, 
where the atomic background energy is chosen to be zero. 
We ignore inelastic collisions - which is a reasonable ap- 
proximation at low density, and let k = Un, where k is 
the bare atom-atom coupling due to s-wave scattering. 

Here a momentum cutoff is implicitly assumed, since 
in renormalizable theories one expects to obtain finite 
results only after the infinities are absorbed through a 
redefinition of bare couplings. To manipulate integrals 
that a priori are divergent, we regularize them by a sim- 
ple cutoff: integrals over k are restricted to |k| < K. 

The homogeneous Hamiltonian, Eq. (J2J), has an exact 
eigenstate in the simplest two-particle sector 00- I n 
momentum space, we expand the field operators 4 f o( x ) 
and ^i(x) in terms of Fourier components a(k) and 6(k), 
respectively, with commutation relations [<z(k), a^(k')] = 
[S(k), St(k')] = S(k - k'). Including a cutoff K, the (un- 
normalized) two-particle eigenstate corresponding to the 



2 



zero center-of-mass momentum is given by 0, Q 



a\0) 



K 



d 3 kG(k) 

(2tt)3/2 



-i W/2 



S+c^S^-k) 



|o), 

(4) 

where N = 2 in the exactly soluble two-particle case, and 
G(k) is the atomic pair correlation function in Fourier 
space. This coherent superposition of a molecule with 
correlated pairs of atoms can be viewed as a dressed 
molecule. More generally, this is also a useful low-density 
variational ansatz for N > 2 particles, where it describes 
a BEC of dressed molecules @ • 

Fermionic or heteronuclear case. Next, we wish to con- 
sider the important case of Fermionic atom pairs (with 
mi = 77i2 = m) in different spin states combining into 
a Bosonic molecule. This is especiallv relevant to the 
studies of ultracold Fermi gases 0, 0, E| in the region of 
resonant superfluidity and BEC-BCS crossover. These 
experiments are notable for the greatly reduced inelas- 
tic loss rate from atom-molecular collisions, due to Pauli 
blocking j^. In this Fermionic case, we only have an 
s-wave coupling between the different fermions, so that 
k = U\2- In addition, the final results of this section can 
be applied to heteronuclear molecules (with either statis- 
tics of the constituent atoms), except that the mass m 
has to be replaced by 2m r , where m r = m,\mil (mx+m?) 
is the reduced mass. 

The Hamiltonian relevant to this case, also has an 
exact eigenstate in the two-particle (N = 2) sector |9j- 
Expanding the field operators \l/i,2(x) in terms of Fourier 
components 6i,2(k), the eigenstate is now given by 



-i N/2 



« T (°)+/ Vw7 b K k ) 6 2(-k) |0) . 
J\k\=0 (27r) J / 2 

(5) 

As before, this is also a useful variational ansatz for the 
TV-particle (N > 2) ground state, where it extends BCS 
theory to include a coherent molecular field. 

Exact eigenvalues. In either the homonuclear or het- 
eronuclear case, the exact energy eigenvalue correspond- 
ing to the two-particle eigenstate (N = 2) is known [1,0. 
Introducing a multiplicity parameter s, where s = 1 for 
the homonuclear case, and s — 2 for the Fermionic or 
heteronuclear case, we find that 



E — En 



Hsx 2 



2ir 2 hra/m 



tqK — tan 1 {tqK) 



(6) 

For real and positive tq, this corresponds to a bound state 
with negative energy, and the resulting binding energy is 
Ei, = —E. The quantity tq is the correlation radius or 
the effective size of the dressed molecule. The right-hand 
side of Eq. © needs to be solved for r , or equivalently 
for the binding energy Eh = h 2 /(mr 2 ,) as a function of 
E m , but in general it has no explicit solution. 



Next, it is useful to re-express the bare Hamiltonian 
parameters in terms of renormalized observable param- 
eters that are invariant at large momentum cutoff. We 
therefore include a nonperturbative renormalization us- 
ing integral equation methods from scattering theory 
1 1 0| . which has some subtle features. In particular, a 
repulsive contact potential with k > has no effect - 
it does not lead to scattering in three-dimensional field 
theory. However, either positive or negative scattering 
lengths can be generated from the same type of attrac- 
tive contact potential with k < 0, depending on the lim- 
iting procedure: if it is carried out with sufficiently deep 
potentials to allow a bound state to form in the atomic 
field channel, then a positive scattering length is possible 
even with an attractive short-range potential. 

The renormalization 01 expresses the bare values as 
k = Tk , x = Txo, and E m = E + shfiTxl/2, in 
terms of the observed or renormalized values kq, xoi 
and Eq. In the Feshbach resonance case, for definitc- 
ness, Eo = Afi(B — Bq). Here, the cutoff K is in- 
cluded through a scaling parameter T = (1 — /3ko) , 
where /? = mK/(2Tr 2 h) 1 kq — Airhabg/m, and a^ g is the 
background s-wave scattering length for the atoms. In 
the homonuclear case, A/j = 2[i\ — fj, m is the magnetic 
moment difference between the atomic and the bound 
molecular channels, and Bq is the magnetic field corre- 
sponding to the resonance, while in the heteronuclear 
case, A[i = ^i + A*2 - A*™- 

We now wish to rewrite Eq. JHJ) in terms of the renor- 
malized constants xo> K o, & n d Eq. After taking the limit 
of large momentum cutoff K, we obtain the following 
simple analytic result: 



Eq = —Eh 



sChx 2 Q VE b 



1 - 2C7ioVE b 



(7) 



where C = m 3 / 2 /(87rft 2 ). Using Eq = Afi(B - Bq), 
Eq. (J2J) can also be rewritten in terms of the magnetic 
fields, so that the resulting binding energy can be di- 
rectly compared with the experimental data. 

In the JILA 85 Rb experiments 0, |2l|, the creation 
of a homonuclear dressed molecular state at a given B 
value is followed by a rapid change in the magnetic field 
which allows an interference fringe to be observed in the 
remaining total number of atoms. The reason for the 
fringe is due to the fact that in a dynamical experiment, 
paired atoms in the dressed molecular wave function can 
interfere constructively or destructively with condensate 
atoms that are not in the ground-state wave function. 
The resulting interference pattern oscillates with a fre- 
quency corresponding to the dressed molecular binding 
energy. 

A graphical solution of Eq. 0, i.e., the binding en- 
ergy E b vs B for the JILA 85 Rb experiment EH], is 
plotted in Fig. ^ together with experimentally observed 
Ramsey fringe frequencies, which are interpreted in the 
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FIG. 1: Binding energy Et versus the magnetic field B for the 
85 Rb experiment The parameter values are taken from the 
subsequent high- precision measurements on the same system 
[2l|: Bo = 155 G and db g = — 443ao, where ao is the Bohr 
radius. In addition, we take xo = 2.84 x 10 -4 m 3//2 /s and 
A/i = — 2.23/is H3], where hb is the Bohr magneton. The 
solid line is our theoretical result, Eq. Q, while the dashed 
line is the result of Eq. J5J . The circles are the experimental 
dataofRef. El. 
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FIG. 2: Binding energy Eb vs B for the 40 K experiment 3, 
where Bo = 224 G and at g = 174ao. In addition, we take 
Xo = 1.12 x 1(T 4 m 3/2 /s and A/i = 1.27/xb 23]. The solid 
line is our theoretical result, while the circles and the error 
bars are the experimental data from Fig. 5 of Ref. Q. 



experiment as a dressed molecular binding energy. The 
agreement between this simple analytic result and the 
experimentally observed binding energy (as well as the 
coupled-channel calculation [10j) is excellent. 

Binding energy measurements were also carried out for 
the case of Fermionic 40 K atoms in two different spin 
states combining into a Bosonic molecule 0. In Fig. |3 
we plot the solution to Eq. J7J), i.e., the binding energy 
Et, vs the magnetic field B, for this experiment, where 
we also see a good agreement between the theory and 
experiment. 

Near-threshold physics. There are common features 
with either Fermionic and Bosonic atoms. All results are 
expressed in terms of the four observable parameters xoi 
kq, A/Lt, -Bo, and are clearly independent of the cutoff, as 



one would expect from a renormalizable theory. There 
are two cases corresponding to different signs of kq: 

(i) Attractive case. If a& 9 < 0, then Eb is a single- 
valued function of B, so there is only one solution branch. 

(ii) Repulsive case. If a& 9 > 0, then Eb is a double- 
valued function of B, so there are two solution branches. 
This case has a bound state in the atomic channel. 

In all cases, the physics near threshold is crucial to un- 
derstanding either type of experiment. For small binding 
energies in the vicinity of the resonance, Eq. {JJ gives a 
quadratic dependence of Eb on Eq (or on the magnetic 
field B): E b ~ Eq 2 / {sChxl) 2 , where E = Afj,(B - B ). 
This is in agreement with the simple resonant scatter- 
ing theory result that Eb = h 2 /ma(B) 2 near the reso- 
nance |l8|, |24|. Here, the effective scattering length is 
a (B) = a bg [1 - AB/{B - B )}, AB is the width of the 
resonance, and the atom-molecule coupling xo c & n be ex- 
pressed via AB as follows: xo — \/8^ab g A^iAB/ (sm). 

In the opposite limit of large Eb, i.e., for magnetic 
fields far away from the resonance, the same equation Q 
gives linear dependence of Eb on Eq (and hence on B) as 
expected, Eb ~ — Eq + shx 2 /(2Ko), including a constant 
energy shift. This linear behavior is not accessible with 
the resonant scattering theory result of Eb = h 2 / 'ma(B) 2 . 

For C\Ko\^/Eb <C 1, i.e., either for small background 
scattering |a& a | or small binding energies Eb near the res- 
onance, we can neglect 2CnQ\fEb in the denominator of 
the second term in Eq. (JJJ) and obtain a quadratic with 
respect to \fEb~. This has the following explicit solution: 



E h 



-En 



4£ n 



sChxl) 2 



(8) 



which (for s = 1) coincides with Eq. (21) of Ref. [l7j . 
This result formally incorporates the above-mentioned 
quadratic dependence of Eb on Eq near the resonance 
where 4Eq/(sCHxq) 2 ^ an< ^ the linear dependence far 
away from the resonance. The quadratic dependence is 
in qualitative agreement with the behavior found from 
our exact result. However, the linear part [Eb — — Eq) 
- while giving the correct slope of the binding energy - 
does not account for the energy shift term sHxq/(2kq) 
due to the renormalization of E m . This leads to a dis- 
crepancy seen in Fig. ^ (dashed line) away from the 
resonance, and is due to the fact that the assumption 
of C\nQ\^fEb -C 1 used to obtain Eq. © is inconsistent 
with the case of large binding energies under considera- 
tion. We note here that it is also possible to obtain the 
exact result from the molecular Green's function method 
of Ref. 0|) ^ the relevant self-energy term is included 
without approximation |25j . 

The relative fraction of the atomic and molecular com- 
ponents in the two-particle eigenstate can be calculated 
using the conserved total number of atomic particles, 
N = Ni + 2N (or N = Ni+N 2 + 2N in the Fermionic 
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case) . At low density, the closed-channel molecular frac- 
tion is (including a factor of 2 to reflect the fact that each 
molecule consists of two atoms): 



2N /N = (1 + 2F/s)~ 



(9) 



Here, F = J c? 3 kG 2 (k)/(27r) 3 , and the correlation func- 
tion G(k) is given by a Lorentzian 



G(k) = G / (1 



(10) 



where k — |k|, Go = sitixoTq / [2h(a,b g — ro)], and we have 
already taken into account renormalization. Taking the 
integral in F we obtain that 



F = 



r 2 



i 3/2 xt 



8?rrg 32irhy/Eb 



1 - 



7 3 / 2 



K 



4irh 2 



'E, 



(11) 



where we have expressed ro = tij \JmEb, so that the final 
result can be analyzed as a function of the magnetic field 
B using the solution to Eq. J7J). 

Combining this result with Eqs. (7J and we find 
that the average fraction of bare molecules in the closed 
channel is typically very small near the resonance. For 
example, for 85 Rb parameters [2lJ it is no higher than 
2N Q /N ~ 0.07, for magnetic fields from B to B ~ 160 G. 
This implies that the structure of the dressed molecules 
and the underlying physics near the resonance is dom- 
inated by the correlated atom pairs rather than by the 
closed channel molecules. 

We can also calculate the atomic pair correlation in 
coordinate space. This is the inverse Fourier transform 
of G(k) given by g(x) = Gq exp(-|x|/r )/(27rr§|x|), for 
|x| > 0. Since ro = h/y/mEb, it is clear that, near thresh- 
old, the bound states are superpositions of molecules with 
pairs of atoms at very long range. Here, one can expect 
modifications |l8j of the binding energy due to mean- 
field many-body interactions of the correlated atoms, 
which have a character similar to Cooper pairs. Such de- 
partures from the predicted binding energies are indeed 
observed [2l| in high-precision Ramsey spectroscopy for 
85 Rb. Similarly, recent collective-mode spectroscopy in 
6 Li has revealed BCS-like behavior near threshold with 
reduced mode frequencies j2|| , quite different to that ex- 
pected for a conventional molecular BEC. 

In summary, a relatively simple field-theoretic model 
for Feshbach coupling has exact solutions for the eigen- 
states in the low-density two-particle sector. It is able to 
accurately predict Feshbach dressed-molecule binding en- 
ergies, and also gives a physical understanding of the type 
of correlated atom-molecular structure produced in these 
experiments. The model has a simple, universal charac- 
ter, and can be used to describe a variety of cases with 
both positive and negative background scattering length. 
Having analytic solutions of the actual eigenstates pro- 
vides an alternative picture that aids in understanding 
these interesting experiments, and is readily usable as a 
starting point to a more complete many-body theory. 
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